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An  Extended  Finite  Element  Method 
for  Contact  in  Multi-Material  Eulerian 

Hydrocodes 

DavM  X. 


ReseaitQbi  supported  by  Awsjf'  Research 

DAADt9-^--t-ISK 


Objective 

Eliminate  the  last  major  limitation  of  multi-material 
Eulerian  FEM  and  FD  formulations  by  the  development 
of  a  general  contact  methodology  for  handling  contact 
with  slip  and  separation. 

Current  mixture  theories  approximate  perfectly  bonded 
materials. 

Bonding  often  masked  by  low  strength  or  material 
failure  at  the  interfaces,  e.g.,  ballistic  penetration. 

Desirable  attributes  for  the  contact  formulation: 

Can  be  implemented  in  existing  codes  without  major  changes. 
Local  algorithm. 

1 )  No  requirement  for  global  topological  information. 

2)  Ease  of  parallelization. 

Allows  for  evolving  topology  —  a  major  attraction  of  the 
Eulerian  formulation. 


Overview  of  Eulerian  FEM 


Mass; 


Momentum: 


pii  =  V  •  CT  ph 

Energy: 


e  =  :  €  -  V  •  (kVT) 


Lagrangian  Finite  Element  Methods 

The  standard  in  solid  mechanics  for  industry. 
Implementation  of  the  Galerkin  method  on  a  computer. 

Based  on  the  weak  form  of  the  eonservation  equations. 

The  eomputational  mesh  deforms  with  the  material. 

Extensive  literature  ranging  from  mathematies  to  applieations. 

Advantages: 

Aeeurate 

Effieient 

Disadvantages: 

Limited  to  problems  with  moderate  deformation. 

Topology  is  fixed. 
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Evilerian  Finite  Element  Methods 

Evolved  from  Finite  Difference  hydrocodes. 
Advantages; 

Handles  arbitrarily  large  deformations. 

New  free  surfaces  evolve  in  a  natural  manner. 

Phase  changes  and  chemical  reactions  are  readily 
incorporated. 

These  strengths  make  the  Eulerian  formulation  attractive 
for  micromechanics  and  other  engineering  applications. 

Disadvantages: 

Computational  cost  (impact  reduced  by  commodity 
computing). 

Diffusion  of  solution  (but  far  superior  to  rezoning 
methods  based  on  interpolation). 

Material  interfaces  are  perfectly  bonded. 


No  longer 
major 
issues. 


Eulerian  Governing  Equations 


Mass: 


Momentum: 


dp 

dt 


+  V  •  {pu) 


H"  ^  '  (p'^  (S)  w)  =  V  •  cr  +  6 


Energy: 

de 


V  •  (ctt) 


=  cr  :  e  —  V  •  (kVT) 


Operator  Splittinj 


Conservation 

Equation 

d(f) 

dt 

+ 

V  •  = 

=  0 

1. 

Lagrangian 

Step 

d4> 

dt 

— 

=  0 

2. 

Eulerian 

Step 

d4> 

dt 

+ 

V  •  {(jyu)  - 

=  0 

Step  1. 

Step  2. 

<t>At 


(j)uAt 


^At 


Ax 


Ax 


(puAt 


TRlfiTHS 
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Explicit  Lagrangian  Step 

Central  difference  integration  in  time 
(standard  practice  for  explicit  FEM). 

ljn+1/2 

^n+1 


=  M-l  {FLterna!  “  / 
=  +  Ata 

=  a;"  + 


Remarks: 

M  is  diagonal  ->  no  linear  equations  to  solve. 
At  is  limited  by  stability. 

Implicit  integration  also  implemented. 
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1  -D  Explicit  Eulenan  Step 


+  —  fj,— 


i+1/2 


i+l/2t+l/2^ 

V- 

original 


i+1/2 


in 


out 


Volume 

transported 

into 

element  / 


Volume 
transported 
out  of 
element  / 


Lagrangian  Stepi  produces  mean 
values  of  solution  variables. 


Donor  Cell  Transport:  uses  mean 
values  but  only  order  accurate. 


Van  Leer  MUSCL:  calculates 
monotonic  slopes  from  mean  data.  2“^  order 
accurate  in  smooth  regions  &  order 
monotonic  at  discontinuities. 
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Two-Dimensional  Eulerian  Step 


Most  transport  algorithms  based  on 
alternating  direction  methods 
(operator  splitting). 

Simple  method  for  comer 
coupling. 

Difficult  to  implement  on 
unstmctured  meshes  (alternating 
spatial  directions). 


Unsplit  methods  are  gaining  in 
popularity. 

More  complicated. 

Greater  cost. 

Better  accuracy  with  careful 
implementation. 


Example  Eulerian  Applications 
Where  Contact  is  Important 


TRITITHS 
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Direct  Numerical  Simulation  of  a 
Metallic-Intermetallic  Composite 

Material  prepared  and  tested  by  Ken  Vecchio. 

Both  projectiles  are  deflected  by  the  material  from  their  original 
trajectories. 

Fundamental  issue:  What  causes  the  deflection? 

Laminate  structure 
Angle  of  impact 


WHA  Penetrator 


MIL  Before  Test 


.308  Penetrator 


leon 


Experimentally  Acquired 
Microstructures 


Generating  realistie  synthetie 
mierostmetures  is  frequently 
ehallenging. 

The  logieally  regular  mesh  of  an 
Eulerian  hydroeode  simplifies 
importing  digitized  images. 

Example:  Metal  matrix 
composite  of  aluminum  and 
alumina  inclusions  subjected  to  a 
compression  test.  Experiments 
performed  by  Prof  K.  Vecchio, 
UCSD. 


Microstructures 


Experiment 


Mesoscale  Simulation 


Analysis  was  performed  using  the  contact  mixture  theory. 
Results  with  simpler  theories  were  unsatisfactory.  While 
there  is  contact  and  separation,  the  relative  slip  between 
surfaces  is  relatively  small. 


X-FEM  for  Contact 


Model  Contact  Problem: 
Taylor  Bar  Impact 


Elastic-plastic  bar  strikes  elastic 
target  at  300  m/s. 

Plane  strain  for  maximum 
mushrooming  of  the  projectile. 

Projectile  properties: 

p  =  1.71  G=.085 
O^  =  2xl0-^ 
A=lxl0“^ 

Target  is  elastic  with  density 
and  elastic  properties  10 

times  higher. 


Mixture  Theories 

Mixture  theories  are  applied  to  elements 
eontaining  more  than  one  material  (“mixed”  or 
“multi-material”  elements). 

Partition  the  strain  increment  between  materials. 
Evaluate  stress  for  each  material. 

Calculate  element  stress  from  the  material  stresses. 

All  interactions  between  materials  occur  in 
mixed  elements  — >  the  mixture  theory  has  the 
role  of  a  contact  algorithm  in  traditional  multi¬ 
material  Eulerian  formulations. 


TRITITHS 


Mean  Strain  Rate  Mixture  Theory 


Timers  .OOE+01 
CycLe=  762 


Bar  Impact  —  EuLerian 


r'l  C'l  m  m 


Contours  show 

original 

coordinates. 
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Zero  Deviatoric  Stress  in  Mixed  Elements 


Timers  .OOE+01 
CycLe=  762 


Bar  Impact  —  EuLerian 


Contours  show 

original 

coordinates. 


Contact  Mixture  Theory 


Timers  .OlE+01 
Cyc[e=  771 


Bar  Impact  —  EuLerian 


r'l  C'l  fn  m 


Contours  show 

original 

coordinates. 


TRITITHS 

UC  GAN  DtESO 


Stress  (MPa) 


What  is  Wrong  With  the  Mixture  Theory? 


Elastic-Plastic  Projectile 


Elastic  Target 


Projectile 


Target 


Distance  (cm) 


slip)  and  a 22  (limited  by  lowest  yield  stress  and  continuous 
across  interface)  are  near  zero  at  the  interfaee. 


GiilS 


Limited  by  the  yield  stress  in  the  projectile. 

Large  in  the  elastic  target. 

Strain  rate  D^j  is  the  same  for  both  the  target  and  projectile  (not  governed 
by  jump  conditions). 


Coupled  Rigid  Body  -  Eulerian 


Contours  show 

original 

coordinates. 

Can  be  improved 
significantly  by 
changing  the 
penalty  stiffness, 
spacing  of  nodes, 
etc. 


New  Eulerian  Contact  Approach:  X-FEM 

Loml  Ewwhmmt  of  tho  Volooity  FioM 


Results  shown  here  are  preliminary. 

Developed  by  Belytschko  et.  al. 

Degrees  of  freedom  g-  added  locally  to  enhance 
accuracy. 

Typical  application  is  fracture.  Enrichment  funetions 
approximate  near  tip  strain  field  and  free  surface  of 
crack. 


Generates  compatible  displacements. 


X-FEM:  Local  Enrichment  of  the 

Velocity  Field 


Nodes  surrounding  all 
the  mixed  elements  are 
enriehed  to  inelude  for 
eaeh  material 
independent: 

Accelerations  &  Forces. 

Velocities. 

Coordinates. 

Masses. 

Mixed  nodes  are  handled 
with  data  structures  and 
algorithms  similar  to 
those  for  mixed 
elements. 


X-FEM:  Local  Enrichment  of  the 
Eulerian  Velocity  Field 

At  the  beginning  of  time  step  n,  the  Lagrangian  coordinates  for  each 
material  coincide  with  the  Eulerian  mesh  coordinates. 

With  an  element,  the  velocity  and  displacement  of  each  material  m 
is  interpolated  independently. 

u”  =  ^Na{x)vI^ 


The  strain  rate  for  each  material  is  calculated  using  its  own  velocity 
field.  ^  ^ 


The  stress  for  each  material  is  updated  using  the  material’s  strain 
rate. 

Unlike  other  X-FEM  applications,  need  to  enforce  contact 
constraints  for  compatible  displacements. 
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Enforcing  the  Contact  Inequality 

Objectives: 

Response  normal  to  the  contact  surface  is  unchanged  from 
standard  Eulerian  formulation. 

Decouple  response  in  tangential  direction. 

Motivate  contact  enforcement  in  the  normal  direction 


based  on  a  one-dimensional  model  problem. 


Enforcing  the  Contact  Inequality 

Multiph  Dmemions 


pTOT 


P=n®n 


F-L  =  [/-P] -F*" 
F^+F-^ 


Separation  Condition 

IFn*o*n>  0  — ^/‘'”  =  F"* 


When  surfaces  separate, 
void  must  be  introduced. 

The  intersections  of  the 
interface  with  the  element 
boundaries  are  calculated  in 
terms  of  the  isoparameteric 
coordinates  (s',?)  at  the 
beginning  of  the  time  step. 

The  volume  generated  by 
the  surfaces  at  the  end  of 
the  time  step  is  calculated, 
and  the  appropriate  amount 
of  void  material  is  added  to 
the  element. 


for  the  Eulerian  Step 


Beginning  of  step 


^ 


Comptiison  of  Lagrangian  and  Bnlorian  Solutions 


Time  =  79.994 

Contours  of  Effective  Rastic  Strain 
max  ipt.  value 
min=0,  at  elem#  1 
max=1 .97096,  at  elem#  1882 


Fringe  Levels 
2.000e+00 

1 .800e+00  I 
1  .GOOe+00  J 

1 .400e+00  _ 
1 .200e+00  _ 
1  .OOOe+00  _ 
8.000e-01  _ 
G.OOOe-01  _ 
4.000e-01 


1 


Time=8.00E-i-01 
Cycle=  764 


Bar  Impact  --  Eulerian 


Eplastic 

2.00 

1.80 

1.60 

1.40 

1.20 

1.00 

0.80 

0.60 

0.40 

0.20 

0.00 


O  Q  O 

o  tn 
d  rn  rn 


Overall  response  is  similar. 

Initial  spatial  resolution  is  similar:  18  x  50~70. 
Spatial  resolution  at  late  times: 

Lagrangian  is  finer  vertically. 

Eulerian  is  finer  horizontally. 

Peak  plastic  strains  are  close. 


X-FEM  at  an  Angle 

Does  the  X-FEM formulation  only  work  when  the  contact 

surface  is  aliened  with  the  mesh? 


Powder  Compaction 

Large  number  of  materials 
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Time=4.50E-01 
Cycle=  1485 


Stainless  Steel  U_p=  1.000 
km/s 
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Ballistic  Penetration 


Steel  projectile. 

340  m/s  initial  velocity. 
1  cm  diameter. 

Copper  target. 

1  cm  thick. 

Mesh  80  X  80. 
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Ballistic  Penetration 


Standard  Mixture  Theory 


New  Contact  Formulation 


\  \ 


Plastic  strain  higher.  \ 

^Greater  penetration 
Projectile  tip  blunted.  depth. 


Conclusions 

Mixture  theories  are  limited  in  their  ability  to  model  contact  and  slip 
due  to  the  limited  information  in  a  single  velocity  field  across  all 
materials. 

Penalty  methods  can  give  good  results  but 

Require  traeer  partieles  for  enforeing  the  eontaet  inequality  eonstraint. 

Suffer  from  the  same  limitations  as  Lagrangian  representations  of  the 
boundary  in  Eulerian  hydroeodes. 

Enriching  the  velocity  field  via  X-FEM  to  give  each  material  an 
independent  velocity  field  shows  great  promise. 

Can  be  retrofitted  into  legacy  codes. 

Uses  only  local  information  for  easy  parallelization. 

Small  memory  burden. 

Eliminates  the  mixture  theory  —  materials  aren’t  artificially 
equilibrated. 

Can  be  generalized  to  include  friction  and  slip  without  separation, 
e.g.,  grain  boundary  slip  in  nanocrystalline  materials. 
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